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ABSTRACT 


The object of this theoretical study is to determine 
the overall magnitude of the effect of an initial stiffener 
web tilt upon the bending stresses induced in a flat plate 
panel subjected to a uniform lateral load. 


A simplified plate geometry was set up, and analyzed 
by use of the theory of small deflections. The theoretical 
solution was programmed for the IBM 7094 computer. Input 
data was selected to represent ship bottom plating, and 
ranged from an aspect ratio (a/b) of 2.0 to 4.0. 


The study revealed that the effect of such an initial 
tilt is significant, and further that the magnitude of the 
effect is essentially independent of the support conditions 
assumed for the plate panel. In addition, a useful design 
parameter, the depth ratio of the stiffener to plating 
thickness (¢/t), was disclosed, and the variation of the 
web tilt effect with this parameter shown. 


Thesis Supervisor: J. Harvey Evans 
Titles Professor of Naval Architecture 
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I. INTRODUCTION 


Background 


"Web Tilt,” as used in this thesis, is defined by 
NAVSHIPS 250-637-3 [1] as "the deviation of the web of a 
stiffening member from a reference plane normal to the base 
line." In recent years, with the increasing importance of 
submarine pressure hull design, the problem of determining 
the magnitude of the effect of web tilt upon a stiffened, 
circular cylinder has come into direct focus. Wenk and 
Kennard [2] and Kennard [3], at the David Taylor Model Basin, 
have studied the effect of such a tilt upon the twisting of 
the stiffening member itself, with the attendant weakening 
of the overall structure. 

This approach to the problem implies that the twisting 
of the stiffening members observed in cases of general in- 
stability failure might have been at least a partial cause 
of the failure rather than simpiy a resultant deformation. 
Despite these theories, there remains an absence of any good, 
experimentally verified theory covering the effect of such 
a tilt upon the stresses obtained in the shell plating it- 
self. In other words, if web tilt did in fact contribute 
to failure, by just how much and in just what way? Asa 
result of this lack of theoretical or experimental informa- 
tion, severe building specifications have been adopted to 
limit the amount of web tilt acceptable in submarine con- 


struction. 
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While the possible existence of a stress magnifi- 
cation which does not lend itself to analysis is of obvious 
concern to the submarine designer, whose structural design 
is a critical art, the possibility that some degree of stress 
magnification exists in shell plating due to web tilt is in 
no way limited to the cylindrical geometry. Virtually all 
ship hulls are made up of stiffened plate panels, both flat 
and curved, and a brief examination of representative plans 
indicates that the orientation of hull stiffeners with 
respect to the hull plating is by no means constant in any 
one ship. If the definition of web tilt given before were 
re-worded such that the angle of tilt was described by the 
angle between the web of the stiffening member and the line 
of action of the normal, or hydrostatic, loading on the 
plate, then the significance of the above statement becomes 
clearer. On some ships very large angles of tilt are formed 
intentionally by using one of the stiffening members near the 
turn of the bilge as a support for some internal part of 
the ship, such as a reduction gear. 

in short, the presence of stiffened plate paneis pos- 
sessed of some degree of web tilt can be assumed on nearly 
every ship. Only the seriousness of this situation, that 
is, the value of any resulting stress magnification, re- 
mains a question. It is to this question that this thesis 


is directed. 
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Statement of the Problem 
Reduced to its simplest form, the "tilt" problem 


arises from the fact that the line of action of the reaction 
force imparted to the stiffening member from the plate under 
conditions of lateral load does not pass directly through 
the flexural center of the stiffener. As discussed by 

Roark [4], this force (P) can be resolved into an equal end 
Parallel force (P') passing through the flexural center (Q) 
of the beam, and a couple (T) ecual to the moment of P and 
Q. As a result of this, an equal and opposite couple (T') 
will act on the plate from the stiffener, causing additional 
stresses and altering the effect of the hydrostatic load 


acting alone. (Figure (1)) 


Ayr‏ رار ےار 





Fig. I Schematic of a tilted stiffener 


Brief consideration of this problem as applied to a 
ship hull reveals that the actual performance of a tilted 


stiffening member under load is a complicated function of 
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many variables. Considering the major variables in turn, 
they are 

l. Accurate description of the boundary conditions 
at each edge of any plate panel. 

Traditionally, ship plate problems have been treated 
as either fully clamped at the edges or else simply sup- 
ported. In an investigation of this sort, however, where 
the action of the stiffening member along one or more of the 
edges must be considered in detail, these approximations ap- 
pear too general. If, as a first approximation of the prob- 
lem, analysis in this thesis is limited to the case of a 
longitudinally framed ship, and the web tilt is considered 
to be applied along one of the longitudinal stiffeners, and, 
in addition, if attention is directec chiefly to the effect 
at or near the center of the plate panei, then the support 
conditions along the transverse edges of the plate panels 
might be reasonably assumed to approximate those of simple 
support. Use of this approximation represents a worst case 
of stress at the mid-panel section. 

Along the longitudinal edges, it can be argued that 
the simple beam deflection of the stiffening member must at 
all times equal the deflection of the plate to which it is 
attached, and also, the rotation of the beam and the plate 
edge must be equal at the point of attachment. These con~ 


ditions describe the conditions of elastic support described 
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by Timoshenko [5] and others, and seems most appropriate for 
this case. The selected boundary conditions will be dis- 
cussed in more detail in the following pages. 

2. Accurate determination of the magnitude of the 
force transmitted from the plate to the stiffener under 
lateral load. 

An exact analysis of this force can only be done for 
a few individual cases, as it depends upon the panel dimen- 
sions, the lateral load, the support conditions, and, in 
some manner, upon the number of adjacent panels between theo- 
retically rigid supports, This last statement can be visu- 
alized if one makes an analogy to a simple beam continuous 
over several supports and subjected to lateral load. Tabu- 
lated solutions for this simple case [6] show that the value 
of the reaction force at a central support will vary by a 
factor of two or more as the number of spans increases from 
one to four. 

34 xact behavior of the tilted stiffener under load. 

Consideration of Figure II will indicate that a 
stiffener with an initial angle of tilt will act in such a 
way as to change the angle when deformed. The amount of 
this change would be different for different points along 
the length of the stiffener. This deformation arises from 
the stress differentials set up in the flange of the member 


as the beam is displaced from its initial position. If the 
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initial angle of tilt were not a constant along the length 
of the stiffener, hut changed from a positive to negative 
value, as might occur in the case of an initial tilt caused 
by variation in welding heat along the length, then it could 
even be argued that the change could cause the final tilt 


angle to decrease. 





Fige II Schematic of stiffener action under load 


4. “he effect of loading in the plane of the plate. 

Clearly, as the stiffener is deflected from its 
initial position by the action of the lateral load, the ef- 
fect of any in-plane load becomes pronounced. In considera- 
tion of ship hulls, such loads are always present, and, as 
in the case above, could act to change the initial angle of 
tilt in any manner, depending upon the sign of the load. 

With all of these considerations in mind, a general 
solution to the web tilt problem is clearly impossible. In 
order to make a first approximation of the magnitude of tilt- 
induced stresses, a simplified model of the situation must 


be defined. 















— —— — 
bl ed A o + as 
— TP αν - wa. a 


kak...‏ — — — ا Un Á...‏ را 
o Qo‏ — ادا س ole ' i | < Hrór Ó  arGw k‏ 


— c cm o im la U, 





-— en AA -« 
١ - - - |) ——- — o -— - -—— 
۱ وك د‎ ٠.٠.26 كين‎ © y «dr ل‎ 
u . “ th μι Ιώ mm a 
u É— ΕΕ 5 O a 





- 0. MENT — UN 





— 


To this end, this thesis will analyze a single plate 
panel by use of the theory of small deflections. This ap- 
proximation has been discussed by Bleich [9] and others, 
who state that the deflection of ship hull plating usually 
does not exceed one-half of the plating thickness, and con- 
cludes that linearized theory may be applied. 

As stated before, two edges of the plate panel will 
be assumed to be simply supported, representing the trans- 
verse bulkhead connections. The other two plate edges, rep- 
resenting the longitudinal stiffening members, will be as- 
sumed to be elastically supported by the stiffening members, 
one of which will be subjected to various angles of initial 
tilt. Since a review of the literature revealed no treat- 
ment of the actual behavior of a skew stiffener under the 
action of a lateral load, the approximation will be made that 
the stiffener deflection will be small enough that the tilt 
angle might be considered a constant under the action of the 
lateral load, i.o., Ø = g' in Fig. Il. 

The angle of initial tilt will be assumed to be a 
constant along the edge of the plate. With this assumption, 
the change brought about by the superposition of an in- 
plane load would be negligible, and will be neglected in 
this paper. Thus the investigation will be limited to the 


effect of the tilt on the plate bending stresses alone. 
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With this simplified approach, some idea of the 
order of magnitude of the tilt-induced stresses ns related 
to those arising from the normal pressure alone will be ob- 
tained. This will not represent a completely general solu- 
tion to the problem in any case, in that no "design equation" 
will be provided. Rather, this approach will represent a 
first, crude look into the problem and, if indeed any prob- 
lem is shown to exist, will point the way to the best 


avenues for further research on the subject. 
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II. THEORY 


Nomenclature 





Fig. III Schematic diagram of simplified plate panel 


a plate length in the «direction (inches) 
b plate length in the y-direction (inches) 
EN coordinate directions 

h, t piate thickness (inches) 

w plate deflection (inches) 

E Young's modulus of elasticity (psi) 


(notes a value of 30,000,000 psi is assumed in this 
work) 

ν Poisson's ratio 
(note: a value of 0.3 is assumed in this work) 

q uniform lateral load on plate (psi) 

ë distance from toe of frame web to the flexural center 


of the stiffener (0Q in rig. I) 


E un 
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g angle of tilt of the frame web 

I moment of inertia of the stiffener about the neutral 
axis of bending (inches”) 

D flexural rigidity of plate (Eh’/12 (1-v?) 
modulus of rigidity of the stiffener 

K a factor analogous to the polar moment of inertia 
of the stiffener section, but dependent on the form 


and dimensions of the cross section 


Solution Technique 


The problem of evaluating the tilt induced stresses 
will be undertaken in two separate phases. First, a solu- 
tion will be obtained for the deflection equation of a flat 
plate panel subjected only to a lateral load. From the de- 
flection equation, stresses in the plate can be computed. 
In addition, the force which the plate exerts on the edge 
stiffeners can also be computed. This force, according to 
the theory of small deflections, is the sum of the shear 
force at the plate edge and the rate of change of the 


twisting moment along the plate edge, viz., 


y | 2 Q Mxy i 
M LE C (1) 


This force will be referred to as the “shear pressure” in 


this thesis. 
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In the second phase of the solution, a similarly 
supported flat plate will be loaded solely by a moment along 
its longitudinal edge. This moment will be equal to the 
product of the shear pressure arising from the lateral load 
and the moment arm arising from the web tilt (c SIN Ø). A 
solution for the plate deflection equation under this loading 
will be obtained. 

If, for purposes of this analysis, the system is 
assumed to remain elastic under action of the loads, the 
principle of superposition may be applied and the separate 
Geflection equations added directly. From this combined 
case, the total stress in the plate panel with a tilted 
stiffener may be obtained, as well as factors relating to 
the percentage of stress magnification at any point in the 


plate panel. 


Fundamental Equations 


The basic equation governing the small deflections 
of laterally loaded plates has been given by Timoshenko [5] 


and others, and is, 


We ὦ = g/D (2a) 
Or, 
4 4 4 
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In the discussion of bending of rectangular plates which 
have two sides simply supported, a solution of the following 


form has been proposed by Nádai [7] and Lévy [8]. 


= Ç y sin EXX (3) 
۱ n 2 
Ls 
m=1 
where, 
ga” ， 
يلا‎ 7 — ay cos BEY + p OOX sim BM ac sım I + 
D 
miry mry | 
IS COSH — (4) 


The series constants in the above equation remain to be de- 
termined by the particular boundary conditions of the case 
in question, and by the equation of the deflection surface 
(2b). 

Since the simplified model selected for examination 
in this thesis has been assumed to be simply supported along 
its transverse edges (x = 0, a), this form of solution is 


considered appropriate for use. 


Boundary Conditions 
As stated before, the solution to the problem is to 


be divided into two parts, each of which having its own par- 
ticular boundary conditions. Certain conditions, however, 


will remain constant throughout the probiem, these are 
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a) The deflection and moment of the plate edge at 


a Simply supported edge must always be zero, viz., 


ΠῚ 


ox 
b) The condition of elastic support requires that 
the load on the stiffening member at the plate edge be 
equal to the force transmitted from the plate to the stif- 


fener, i.e., the shear pressure, viz., 


90) 


EI EJ Μα at y = + b/2 (6a) 
x 
or, 
e 3 y 
EI oz D. CARS 6 (29v) —T laty 28 b/2 (6b) 
3x” ay? 32 dy 


Considering other conditions which must exist along 
the plate edges, it can easily be seen that the rotation of 
the stiffener at its point of attachment to the plate must 
equal the rotation of the plate at the same point. Stated 


in mathematical terms, this would be 


33 ۴ 2 ὦ δῷ 7 
OE AA AS AAA | ناه‎ y = + b/2 (7) 
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In most actual cases, however, plating stiffeners 
are T-sections, and it can be shown that the value of the 
torsional stiffness constant (K) is only a fraction of the 
value of the poiar moment of inertia (J) of these sections 
[4]. Hence, the resistance of such a structural element 
to a rotation about its point of attachment is quite small. 
This is obviously not a good approximation at or near the 
ends of the stiffener, where some form of connection is 
made to a transverse member, but near the center of the 
panel, where the bending stresses resulting from boundary 
condition (5) would be the greatest, the assumption seems 
reasonable. 

Considering this, one further assumption which might 
be made concerning the longitudinal plate edge would be 
that the stiffener, near the center of its span, offers no 


resistance to rotation about its point of attachment, or, 


VE e 
$i + ANO at y = + b/2 (8) 
8 dy ax? ^7 


It is just this assumption which has been made by Timoshenko 
[5] in his treatment of a somewhat similar case, There, the 
value of the stiffness constant (K) was equal to the polar 
moment (J), a fact which serves to affirm the assumption 


made in (8). 
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Another approach to the condition existing along the 
+ b/2 edges may be defined if one considers an ideal plate 


deflection pattern as shown in Fig. IV. 





Pig. IV. Ideal deflection pattern of an infinite series of 


similar plate panels under lateral load 


While this deflection pattern is obviously an ideal situa- 
tion in which each plate panel affects its neighbor in such 
a way that all rotation of the piate panel at the stiffener 
is removed, viz., 


Au 


iO T. at y = + b/2 (9) 


it would appear reasonable to assume that this situation is 


no more ideal than that which assumes no residual bending 


moment at the stiffeners (8). Actually, it seems most likely 


that the actual situation lies somewhere between the two 
extremes. Accordingly, both of these approximations will 
be made in the solution of the plate under lateral load 


alone. 
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For purposes of identification, the solution using 
boundary condition (8) will be referreá to as the "simple 
support” approximation to the lateral load case, and that 
using boundary condition (9) as the "clamped" approximation 
to the lateral load case. This terminology will be retained 
throughout. 

The second phase of the solution, that dealing with 
the plate panel loaded only by a moment applied along a 
longitudinal edge, has been discussed by Timoshenko [5] for 
somewhat different boundary conditions. A completely 
general edge moment loading may be described by the fol- 


lowing two conditions: 


- a o وه‎ 7 
οι πο τν πο = f, (x) at y = + b/2 (10) 
dy IX 
2 2 _ 
-οἱ581ν 55 περα) at y = - b/2 (11) 
- ey a 7 N 


Using these conditions, together with (5) and (6), the com- 
plete solution for the plate deflection equation under moment 
load alone may be obtained in general form. Specialization 
of this solution, by setting the moment at one of the longi- 
tudinal edges equal to zero and the other equal to that 
transmitted from the tilted stiffener, will lead to the 
solution sought in this phase. This method of specialization 


is further described in appendix B. 
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IZI. PROCEDURE 


Lateral Load Case 


l. Simple support approximation 


The boundary conditions appropriate to this case are 














2 
v = Û, — 0 at x = 0, a (5) 
3x^ 
A 6. 2 3 نی‎ 
EI Bap S84 (zy) شك‎ aty=+th2 (6) 
ox “dy ox dy ^ 
P E 
ay ax^ 7 


The solution of equation (2a) for these conditions has been 
carried out by Timoshenko [5] and, with a slight change in 


notation, is reproduced below. 


4 cn 

_ 4ga ` <j l mT 

ον πα ` s 1A cosi E 
T D m 
8-21 325 ٠ ه‎ 
a 
+B U sm X | gy BE (12) 
m a a j a 


in which, 


v (1+v) STNH Av (Lev) A, Cosi Ag 70, (2COSH AQ uS EH An ) 


2 2 
(3+v) (1-v)SINH À „COSH ey ) A, +20, COSH A 





toT 


m 
(13) 


v (1-v)SINH سا‎ COSH 5 





B = 
m 2 λ 


(3+v) (1-v)SINH À COSH An (1-ν΄}λ +20, COSH £ 
(14) 


a 


æ) Qe 


= mrb 9 
* ` 2 (15) 
.EL mm 


with this solution for the plate deflection surface, 


surface stresses in the plate may be found from the following 


relationships: 
6M 2 2 
am. = — = سے‎ E Γ ο : + y Ss» J 1 (17) 
= h^ - àx oy” 7 
6M Ta 2f. 
= -3 P... (18) 
Y h h dy ox 


Upon differentiation of (12), and substitution, the 
equations for the stresses due to lateral load under the 


simple support approximation are 


_ 22و24‎ ς i j 


9 my 
Do e ~ 3 ‚tra, (1 v)COSH = 


m x BY | n my Dry ` Dnm 
2v B „COSH : 5 BAL v) š SINH - | SIN = 


(19) 
B - 
σ = 24da | i. | vt (v-1)A COSH τν 
ys hiro L. m b. m a 
πε, ως 


on miry = mmTy mmy ^ y TE 
2B, COSH z + (y 1)B, E SINH 2 J SI = 


(20) 
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2. Ciamped approximation 


The boundary conditions appropriate to this case are 


. = 0, TE 0 at x = O, a (5) 
Ox 
3 E > 
prep) 24 (2-) so | at y = + b/2 (6) 
ðX dy Ox ΟΥ 
AW 
< 一 二 0 at y = + b/2 (2) 
oy 


A detailed solution for this case is presented in Appendix A. 


The solution for the deflection surface of the plate 


subject to the approximation that there is no rotation of 


the plate panel at the stiffeners under purely lateral load 


is 
4 4 c 14 rà mTwv ra Ta Ty 
yg = =a — lsa σοδη 553: +p A srg BX gry BA 
Tp m? L m a m a a 11 a 
۲۳۱21 3 هم‎ 
(12) 
where 
"e ~a (SINH batis COSH A) (a 
m 2 
Anim m INH À „COSH À mt Z SINH iu 
q SINH A 
Br = e ` (22) 
Cry Ay FO, SINH Ay COSH A) +2SINH À 
with 


the terms An and αι as defined above by (15) ana (16). 
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The solution is thus of the same form for both ap- 
proximations. Only the series constants, A, and B,» change 
with different boundary conditions. With these same modi- 
fications then, the equations derived for the flat plate 


stresses are identical for both support approximations. 


Tilt Load Case 





Boundary conditions for this case are 


32 y 
۱ ٩ ۷۶۰ و‎ κ. at x = Q, a (5) 
ox 
at, n 33, E | 
85 >.) f o Re Ἡ at y=+b/2 (6) 
Ox 7" 8y ax ΟΥ 
Ou — Ἴ 
-D ١ — ۷ € = £4 (x) at y = + b/2 (10) 
- dy ax - 
1 XP 32, ۴ 
πι iru SE KE r. (z) at y = - b/2 (11) 
- dy doe = 


A detailed solution for this case is presented in Appendix B. 
If the arbitrary moments applied to the plate edge 


are assumed to be of the general form 


f (x) w ` En SIN — (23) 


then the solution to this case can be given as follows 
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2 «= AE 
y=- ف‎ a! 2 cosı MEY 4 MEY ατνη PTY 
2 ۳۹ 2 LR a aR a 
TD m 1 1 
m=1 
Y m Ty miry mTy 4 TX 
> ' — as 
R, sINH z + aR, COSH 2 1 SIN š (24) 


where the constants of the series terms are given by 


(1*v-a, À) SINH Am (Iy) An COSH A 
— —— — — (25) 
(1--) SYNH سا‎ COSH Am 


l+v=g_ 1 _)C05H A -(1l-=v)A_ SINH A 
(1-v)COSH 入 +C SINH A 
m m τὰ 


r κ 
ز 2 (12۷) و‎ COSH dat tm (17V) SINH hin (27) 


2 
RS 
ii 


— 


R, = | y(lev)+2 


2 ] J SINH Ami. (17y) COSH A (28) 


with the terms λα and a, as defined before, (15), (16). 

In order to determine the value of the series con- 
stant En in (24), the nature of the loading moment must be 
considered. As stated earlier, the moment applied to the 
plate from the stiffener will be the product of the shear 
pressure at the edge of the plate times the moment arm arising 
from the angle of tilt, («c SIN p). 

The shear pressure at the plate edge, by the theory 


of small deflections, is 
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um _ AMxy 
Vz Q, ay (29) 
or 
3 3 
V = =D Í 3 5 * (2-v) 2 1 (30) 
Y - dy ax^ay 7 


The reaction force applied to the plate is then 





3 3 
== Y = D l چ‎ + (2=y) ο | (31) 
y - ay ax“ ay 
y=b/2 y=b /2 


Differentiating the solution to the lateral load case 


(12) and substituting that value into (31) gives 


4qa c bee ue 
Vip WS i — OO) SINH 5 + (v+1)B, SINH An 
m 
m= ۳ 
1 MIX 
v (v-1)B, A COSH ^n SIN 3 (32) 


This solution has the same form for both lateral 
load stress approximations, only the values of A, and B, 
change, as has been shown. 


Since the moment applied to the plate is of the form 


f(x) = 7V (572) e SIN Y (33) 


and has been assumed to be 
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£ ۳ MIR 
£(x) , E, SIN = (23) 


the final fozm of the solution for the series constant E, is 


4qa 
En = 2 2L (v-1)A. SINH Aja F(V+1)B, SINH κ 
m 
+ (v-1)B N COSH Ami c SIN Z (34) 


With this, the solution for the deflection surface 
of the plate panel subjected to both lateral and tilt loadings 
separately is complete. 

It should be noted at this point that when the two 
deflections are combined by the principle of superposition, 
a different form of the equation for the shear pressure at 
the plate edge would arise. This is to say that the loading 
by the tilt moment would give rise to its own shear pressure 
and a separate moment load on the plate. This increment in 
the shear pressure has been analyzed and shown to be second 
order. The analysis is shown in detail in Appendix C. For 
the purpose of this investigation, the second order terms 
are disregarded. 

The stresses induced in the plate subjected to web 
tilt loading may be evaluated in the same manner as before, 


by use of equations (17) and (18). 
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Differentiating (24), and substituting into (17) 


and (18) gives 


3 "o_o (n(v-1)+2v) ۱ _ 
Dec c! — cosa BEY , OD sy TY 
2 h 一 1 a aR, a 
m=l.316.. 


lv (Y-1)+2v) . nry , mmy(v-1) nity ` MTT 
gp AOD ey) STNH BEY + ππγίν-}) cosy 2 emy, E 
R a aR, a 8 a 


2 
(35) 
en 
3 ۲ (p(1-v)+2) Ἴ = à 
OS = 4 Eqs ————— cosg EXY +, Mgr ΞΤΖ 
y a m R a aR a 
2 رد وه‎ 1 1 


1-v)+2 4 
ΟΞ ora € + utm COSH r ¡SIN - 
2 


R, a 


(36) 


According the the principie of superposition, the 
stresses obtained for the two separate cases may be added 
directly to give the total stress in the plate panel sub- 
jected to both the lateral load and the tilt load arising 
from a single tilted frame. 


The total plate stresses are 


J = g + G (37) 


g ع‎ 5 + y (38) 
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and the stress magnification factor arising from the tilted 
frame may be readily obtained as well. 


P4 — P4 — (39) 


F =, ( 40 
9 y 7 %y, (40) 


The procedure shown thus far will determine the 
magnitude of the effect of web tilt for each of the two ap- 
proximations of plate restraint which have been proposed. 

A computer program has been devised which will allow the 
infinite series terms involved in the solution to be evalu- 
ated out to vanishingly small terms. This program is pre- 


sented in Appendix D. 


Selection of Input Data 








Having obtained an analytical solution for the effect 
of web tilt loading on plate bending stresses, some numeri- 
cal examples must be obtained for examination. 

In order to obtain a set of representative input 
data, the design path described by St. Denis [10] was fol- 
lowed through a single cycle. Three values of ship length, 
transverse frame spacing, and longitudinal frame spacing 
were selected as typical values, and a design made for each 


of the twenty-seven possible combinations. The assumed 
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initial values are shown in Table I. 


Transverse 
Spacing 


72 


Longitudinal 
Spacing 


24 





84 


96 






Table I. Initial Input Data 


The value of the primary hull girder stress was as-~ 


sumed to be given by the relationship 


2240°/L (41)‏ = رن 
then, following the procedure outlined in [10], stress‏ 
Schedules were set up for each case to give allowable values‏ 
of the plate bending stress (54). With this, the plate‏ 
thickness required was obtained by using the maximum value‏ 


given by the equations shown below. 


i. plates in pure bending 


ya /0.6859 5H (12) 
A ue 
2 


3 





— 


2. stability of plates in compression 


2 
y Leds da 
t= — — (43) 


AE 


5 [σ᾽ 


The plate thus selected in each case was then checked 
to insure that a factor of safety of 1.25 was provided on 
the ultimate strength after buckling. 

With the plating thus selected, the scantlings of 
the longitudinal stiffeners were selected by setting a limiting 


value of the slenderness ratio (L/k) for each span, viz., 
(L/k) « 30 (44) 


The bending moment in each stiffener was approximated by 


means of the relationship 


2 
۵ HbL 
ue EV πμ (45) 
12 E 


where, in all cases in this thesis, the draft of the ship 
was held constant at tnirtv feet. 
The required section modulus of the combined beam 


and plate was obtained by the relationship 


£ (required) = كل‎ (46) 
C G 
2 
where 
2 (a, / 1.25) - 9, (47) 
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Actual selection of the stiffener sections was ac- 
complished by the use of references [6] and [11]. Medium 
steel construction was assumed throughout. 

A complete tabulation of the results of this first 
design spiral is given in Table II. While it is acknowledged 
that these scantlings are not likely to be appropriate for 
specific ships of the assumed initial dimensions due to con- 
siderations of overall section modulus, it is felt that the 
values obtained in this limited analysis represent a minimum 
design condition, and are thus appropriate for use in the 
analysis for the effect of web tilt. 

In order to provide a more orderly form for data pre- 
sentation, the input data shown in Table II have been pre- 
sented in non-dimensional form in Table III. 

Since the most critical factor affecting the magni- 
tude of the tilt moment can easily be seen to be the Cistance 
from the toe to the flexural center (¢), the input data, 

(a/b) ratios, are shown plotted against the ratio ¢/t)in 


Figure V. 
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300 
300 
300 
300 
300 
300 
300 
300 
300 
500 
500 
500 
500 
500 
500 
500 
500 
500 
700 
700 
700 
700 
700 
700 
700 
700 
700 


72 
72 
72 
84 
84 
84 
96 
96 
26 
72 
72 
72 
84 
84 
84 
96 
96 
96 
72 
72 
72 
34 
84 
84 
96 
96 
96 


24 
30 


24 
30 
36 


30 
36 
24 
30 
36 
24 
30 
36 


30 
36 
24 
30 
36 
24 
39 
36 
24 
30 
36 


Table 





h 


0.375 
0.5 
0.5625 
0.375 
0.5 
075925 
0.375 
0.5 
0.5625 
0.4375 
0.5625 
0.63875 
0.4375 
0.5625 
0.6375 
0.4375 
0.5625 
0.6975 
0.5 
0.625 
0.75 
0.5 
0025 
0.75 
pis 
0.625 
0.75 


ll. 


78.0 
11767 
15127 
107.6 
187.6 
198.7 
158.1 
228.6 
324.8 

90.1 
102,9 
216.2 
13323 
198.7 
296.4 
174.4 
268.2 
32171 
125.8 
176.3 
249.9 
187.6 
25058 
342.0 
228.6 
342.5 
441.6 


Table 





11.4 
1553 
21.1 
19.0 
30.7 
30.7 
30.7 
37.8 
42.2 
13.0 
21.1 
30.7 
VENUE 
1657 
42.2 
30.7 
42.2 
59.6 
1920 
2345 
A352 
30.7 
37,8 
435.6 
37.8 
42.2 
64.8 


5.78 
5.85 
6.72 
0.74 
7.72 
1,72 
7.72 
7.78 
8.72 
5.82 
6.77 
,ω ὁ 
6.77 
7.72 
7.78 
7.72 
7.78 
3.75 
6.74 
6.8 
7.75 
7.72 
7.75 
7.61 
7.78 
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INPUT DATA 


NOMENCLATURE 


long side dimension. 
short side dimension. 
thickness of the plate 


Ios moment of inertia of the stiffeners. 


es 


2.0 


2.33 
2.4 
2.66 
2.8 
3.0 
5972 
3495 


4.0 


distance from the flexural center of the stiffener to 
the plate edge. 
length of the ship. 





10.65 12.0  1.90*10 7 128  12.1*1077 72" 
12.40 27 1.86 149 11.8 84 
12.40 11.7 2.55 144 193 72 
11.00 15.5 1.30 170 M 96 
10.90 15.4 2.24 165 13.4 84 
12.45 15.4 2.46 193 15.8 72 
12.35 15.6 2.67 192 16.0 96 
12.45 18.0 2.59 224 16.5 84 
12.45 20.6 2.79 256 17.8 96 
Table III. Non-dimensional Input Data 
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IV. RESULTS 


The computer programs, Appendix D, were written to 
evaluate each of the infinite series presented as solutions 
for the plate stresses in Chapter III out to its tenth term. 
This extent was considered necessary to insure accuracy in 
the computation of the plate stresses, since the series 
representing the moments in the plate converge less rapidly 
than do those for deflection (i.e., as 1/m? instead of as 
1/m?). 

Figures VI and VII are plots of the results obtained 
for the stresses, in the y direction only, using the simple 
support approximation. Figures VIII and IX show the same 
results obtained by use of the clamped approximation. These 
figures show only the results for case l. As the general 
shape of the resultant curve was identical for all cases, 
differing only in magnitude, presentation of each case in 
this manner was considered superfluous. 

Figures X through XV present a comparison of the 
actual tilt stresses obtained for each of the lateral load 
support approximations used. Again, only selected cases 
have been presented, but these are representative of the 
form obtained in each case examined. 

Figures XVI through XXIII are plots of the results 
of all cases examined, both for the simple support and 
clamped approximations. These plots present the tilt stress 
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obtained, shown as a percentage of the lateral load stress, 
versus the characteristic plate parameters in a non- 
dimensional form, i.e., the depth ratio (c/t) and the aspect 
ratio (a/b). These are shown for three different angles 

of initial web tilt: three, five, and seven degrees. 

In all cases shown in Figures XVI through XXIII, 
data from the simple support approximation has been pre- 
sented at the center of the plate panel (x=a/2, y=0) and 
data from the clamped approximation has been shown at the 
edge of the plate panel (x=a/2, y=b/2). 

Tabulations of the data used to plot the figures 


are presented concurrently in Tables IV through XIII. 
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Figure VI. Sample stress pattern, 
s, stress along x-A/2 axis, Case l. 


Simple support approximation. 
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SIMPLY SUPPORTED CASE (f = 7°) 





These values remain constant for all angles of tilt 


L: 300' 
eft “y no tilt* ὃν tilted % increase 
1 11.7 31,863.3psi. 33,278.2psi. 4,44 
| 12.05 33,115.8 34 ,459 .0 47 
10 35,865.9 37,605.3 4.35 
15.40 33,585.9 35,618.7 6.05 
15.40 38,846.9 41,070.8 5.72 
15.50 37,671.6 39,799.3 5.65 
15.6 34,572.3 36,718.3 6.20 
18.0 39,878.92 42,579.5 6.77 
20.6 40,297.4 43,453.3 7.83 
L: 500' 

DA 7y no tilt ὃν tilted % increase 
IB 2 22,134.7psi. 25,145 563. 4,60 
11.3 23,975.5 25,145.5 4.88 
12,0 25,171.8 26,465.2 SA 
12,7 25,194.3 26,621 .0 5.65 
13.3 28,513.6 30,154.6 5.76 
13.7 26,507.5 30,154.6 6.04 
15.5 29,280.3 31,270.8 6.10 
17.60 29,669.4 31,985.7 7.81 


and are not repeated in the subsequent tables 


Table VI. 


Lateral and tilt stresses 


vs. depth ratio (s/t), 7? tilt angle. 


All cases, simple support approximation. 





L: 700" 
eft 


10.30 
10.30 
10.90 
11.70 
12.45 
13.50 
14.00 
15.40 
15.60 


Sy no tilt 


18,571.5psi. 


20,116.3 
20,364.4 
21,153.1 
21,462.0 
21,829.2 
22,090.3 
22,416.5 


22,708.6 


Table VI. 


— 


Sy tilted 


19,429.Spsi. 
21,076.5 
21,414.4 
22,354.3 
22,769.0 
23,283.9 
23,624.4 
24,154.2 


24,492.0 


% increase 


4.62 
4.17 
5.15 
5.68 
6.09 
6.67 
6.94 
7.75 


7.86 


Lateral and tilt stresses 


vS. depth ratio (c/t), 7? tilt angle. 
All cases, simple support approximation. 








αμ : 


ELIT 

-- -—— y 
Ἱ د‎ 
- uS 


— -——- E سس‎ 


Je tie” و‎ 
=> - — > 








Figure XVI. Stress magnification vs. 
depth ratio (e/t), pe 
Simple support approximation. 
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SIMPLY SUPPORTED CASE 
L: 300' 
G σ 
y tilted, | y tilted, | 

eft (#=3 `) % increase (055) % increase 
11.70 32,47171.0 1.90 32,875.2 3.20 
12.05 23,093,.0 T. 7 34,075.0 2.90 
13.70 36,613.0 2.10 37,110.5 3.30 
15.40 33, 802.0 2.50 40,437.0 4.10 
15.40 34,459.0 3.45 35,040.0 4.33 
25.50 33,585.0 2.43 39,193.2 4.04 
15.60 35,494.4 2.66 36,107.3 4,44 
19.00 ع2 0 م نا‎ 41,810.0 4.34 
20.60 41,652.0 3.36 42,554.4 5.60 
Ls: 500' 

^Y tilted, ὃν tilted, 

s/t (Y=3”) % increase (@=5 `) % increase 
11.20 22 072.0 1098 22,864.0 3.30 
11.30 24,478.4 210 24,813.0 3.50 
12.00 25,727.2 cere 26,096.8 3.70 
12.70 25,807 .2 2.43 26,214.7 4.05 
13.39 29,218.0 2047 29,687.0 4.12 
13.70 27 32953 2.60 27 و 653و‎ 3 4.32 
13.85 28,027.0 2870 28,509.7 4.43 
15.50 244135 1 2.92 30,703.8 4.86 
17.60 30,064.0 3.35 31432999 5.58 


Table VII. Lateral and tilt stresses vs. 
depth ratio («/t), 3? and s? tilt angle. 
Simply supported approximation. 








L: 700' 
Sy tilted, | Ty tilted, 
eft (@=3 `) % increase (@#=5 `) % increase 

10.30 18,490.0 1.98 19, 185.3 3.30 

10.30 20,528.8 2.05 207803.1 3.421 

10.90 20959125 3 2-4 21113 2. 95 

11.70 21,663.9 2.84 200012. 4.06 

12.45 21011 2 2.61 22,396.9 4.35 

13,50 2244539 2.86 22,869.0 4.77 
| 14.00 22, 139,1 24,98 23,187. 5 4.97 
15۰40 23,162.7 3.33 23,659.2 5.54 
15.60 23,475.6 3.38 23,986.0 5.63 


Table VII. Lateral and tilt stresses vs. 
depth ratio (e/t), 3? and 5? tilt angle. 
Simply supported approximation, 
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Figure XVII. Stress magnification vs. 
depth ratio (e/t), 3° and 5° tilt angle. 


Simple support approximation. 
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SIMPLY SUPPORTED CASE (#=7°) 


L: 300' 
A/B y no tilt “y tilted % increase 
2.00 33,115.8psi. 34 ,459.0psi. 4.07 
2.33 35,865.9 37,606.3 4.85 
2.40 31,863.3 33,278.2 4.44 
2.66 37,671.6 39,799.3 5.65 
2.80 33,585.9 35,618.7 6.05 
3.00 38 , 846 .9 4128 ΕΕ 
| 3.20 34,572.3 36,718.3 6.20 
| 3.50 9 42,579.5 6.77 
4.00 40,297.4 43,453.3 7.83 
Ls 500° 
A/B ^y no tilt ^y tilted % increase 
2.00 22,134.7psi. 23,155.0psi. 4.60 
2.33 23,975.5 25,145.5 4.88 
2.40 25,171.8 26,465.2 5.14 
2.66 25,194.3 26,621.0 5.66 
2.80 26,507.5 28,110.4 6.04 
3.00 28,513.6 30,154.6 5.76 
3.20 27,300.0 28,993.0 6.20 
3.50 29,280.3 31,270.8 6.10 
4.00 29,669.4 31,985.7 7.81 


Table VIII. Lateral and tilt stresses vs. 
aspect ratio (a/b), 7° tilt angle. 
Simple support approximation. 
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Le 700' 
A/B Sy. no- tilt “y tilted % increase 
| 2.00 18,571.5psi. 19,429.8psi. 4.62 
| 2.89 20,110.56 2170763 4.77 
| 2.40 20,364.4 21,414.4 5.15 
| 2.66 21,153.1 22,354.3 5.68 
2.80 21,462.0 22,769.0 6.09 
3.00 21,829.2 23,283.9 6.67 
3.20 22,090.3 23,624.4 6.94 
3.50 22,416.5 24,154.2 7.75 
4.00 22,708.6 24,492.0 . 7.86 


Table VIII. Lateral and tilt stresses vs. 
aspect ratio (a/b), 7? tilt angle. 
Simple support approximation. 





πα. 


E 








0 | 


uu 


κε 





-564 


Figure XVIII. Stress magnification vs. 
aspect ratio (a/b), 7° tilt angle. 
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SIMPLY SUPPORTED_CASE 





L:300' 
“y tilted “y tilted 
A/B (g-39) % increase (9=5°) % increase 
2.00 33,693.0psi. L5 34,075.0psi. 2.90 
2.33 36,613.0 2.10 37,110.5 3.50 
2.40 32,471.0 1.90 32,875.0 3.20 
2.66 38,585.0 2.43 39,193.2 4.04 
2.80 39,802.2 2.50 40,437.0 4.10 
3.00 34,459.0 3.45 35,040.0 4.33 
3.20 35,494.4 2.66 36,107.3 4,44 
3.50 41,037.0 2.90 41,810.0 4.84 
4.00 41,652.0 3.36 42,554 .4 5.60 
28500 
Sy tilted Sy tilted 
A/B (0-35) % increase (9259) 9; increase 
2.00 22,574.0psi. 1.98 22,864.0psi. 3.30 
2.33 24,478.4 2019 24,813.0 3.50 
2.40 25 ,727.2 Bou 26,096.8 3.70 
2.66 25,307.2 2.43 26,214.7 4.05 
2.80 27 295.8 2.60 27,653.8 4.32 
3.00 29,218.0 2.47 29,687.0 4.12 
3.20 28,027.0 2.70 28,509.7 4.43 
3.50 30,135.1 2.92 30,703.8 4.86 
4.00 30,664.0 3.35 31,325,3 5.58 
Table IX. Lateral and tilt vs. 


aspect ratio (a/b), 3 


0 


and 5° tilt angle. 


Simple support approximation. 





L:700' 
“y suas | ὃν و لد‎ 

A/B (@=3`) % increase (2=5 % increase 
2.00 18,490.0psi. 1.98 19,185.3psi. 3.30 
2.33 20,528.8 2.05 20,803.1 3.41 
2.40 20,815.3 2.21 21,158.0 3.35 
2.66 21 9 2.44 225018 .1 4.06 
2.80 2 202352 2.61 22,396.9 4.35 
3.00 22,453 9 2 6 22 و‎ 869 . 0 4.17 
3220 22,769. 1 2.98 23920145 4.97 
3.50 23, 108% 7 3.33 23,059.2 5.54 
4.00 23,475.7 3.38 23,986.0 5.63 


Table IX. Lateral and tilt vs. 
aspect ratio (a/b), 3° ana 5° tilt angle. 
Simple support approximation. 
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Figure XIX. Stress magnification vs. 
aspect ratio (a/b), 3° and 5° tilt angle. 


Simple support approximation 
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CLAMPED EDGES CASE (f=7°) 


L:300' 

S 5 f % increase 

eft y no tile y tilted (decrease) 
11.70 23,247.0psi. 20,056.0psi. 13.30 
12.05 26,542.0 23,006.0 13.30 
13.70 26 , 10 22, 0 15.00 
15.40 26,680.0 22 2 120 16.75 
15.40 23,499.90 19,298.0 17.90 
15.50 26 0 220 16.90 
15.60 285423. 29 19,198.0 18.00 
18.00 26,629.0 21,559.0 19.00 
20.60 26,775.0 20,817.0 22.20 

L:500' 

% increase 

e/t “y no tilt 3 tilted (decrease) 
11.20 17,590.0psi. 14,900.0psi. 15.30 
11.30 17,657.0 14,952.0 15.30 
12.00 18,344.0 15,427.0 15.90 
12.70 17,699.0 14,679.0 17.20 
13430 19,401.0 16,102.0 17.00 
13.70 18,366.0 15,053.0 18.00 
13.85 18,374.0 15,040.0 18.20 
SO 19,492.0 15,656,0 19.70 
17.60 19,544.0 5,170.0 22.40 
Table X. Lateral and tilt stresses vs. 


depth ratio (c/t), 7° tilt angle. 


Clamped approximation. 
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Table X. 


L:700' 
% increase 
e/t oy no tilt oy tilted (decrease) 
10.30 17,657.0psi. 14,962.0psi. 15.30 
10.30 14,667.0 12,459.0 15.60 
10.90 14,713.0 12, 346-0 16.10 
11.70 14,761.0 12,219.0 17.20 
12.45 14,814,0 12,112.0 18.30 
13.50 14,834.0 11,918.0 19:275 
14.00 14,700.0 11,678.0 20.60 
15.40 149 11,564.0 22.50 
15.60 14,384.0 11,510.0 25.10 


Lateral and tilt stresses vs. 


depth ratio (e/t), 7° tilt angle. 


Clamped approximation. 
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Figure XX. Stress magnification vs. 
depth ratio (c/t), r tilt angle. 


Clamped approrimation. 
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CLAMPED EDGES CASE 


Ty tilted “y tilged 
eft (g=3~) % increase (d=5") % increase 
11.70 21,877.0psi. 5.90 20,965.5psi. 9.80 
12.05 25,023.5 5.75 24,013.2 9.45 
13.70 24,946.06 6.40 23,802.5 10.75 
15.40 24,762.0 7.20 23,485.0 11.95 
15.40 21,695.0 7.06 20,494.0 12.80 
15.50 24,735.0 7.30 23,449,0 12.10 
15.60 21,607.9 7.74 20,401.1 12.70 
18.00 24,528.3 7.90 23,041.6 13.50 
20.60 24,216.8 9.58 22,514.0 15.90 
1500 
“y tilted “y tilted 
eft (fz3") % increase (gZ=5`) % increase 

11.20 16,435.0psi. 6.60 15,666.0psi. 10.90 
11030 16,499.0 6.60 15,729.0 10.90 
12.00 17 ,091.4 6.85 16,258.0 11.40 
12.70 16 ,402.0 7.30 15,539.0 12.20 
13.30 17,984.0 7.32 17,042.0 12.20 
13.70 16,943.2 745 15,996.5 12.90 
13.85 16,942.0 71.78 15 ,989..9 13.00 
15.50 17,844.0 3.48 16,748.0 14.10 
19560 17,665.0 5323 16,415.6 16.00 


Table XI. Lateral and tilt stresses vs. 
depth ratio (e/t), ΒΡ and 5° tilt angle. 
Clamped approximation. 
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L:700' 
G سس‎ f Y tilted 
eft (053) % increase (0-5) % increase 
10.30 13,663.0psi. 7.10 13,016.0psi. 11.28 
10.30 13,663.0 2.10 13,016.0 11.25 
10.90 13,696.7 6.95 13,020.3 11.50 
11.70 13,669.4 7.50 12,943.0 12.30 
12.45 0 7.84 12,881.0 13.05 
13.50 13,586.8 8.43 10D I.I 14.05 
14.00 13,402.0 8.35 12,539.0 14.70 
15.40 13,474.4 9.67 12,517.0 16.10 
15.60 13,434.4 9.75 12,471.0 16.20 
Table XI. Lateral and tilt stresses vs. 


depth ratio (e/t), g^ 


and 


2 


ES 


tilt angle. 


Clamped approximation. 
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Figure XXI. Stress magnification vs. 


0 


depth ratio (c/t), 3° and 5° tilt angle. 


Clamped approximation. 
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CLAMPED EDGES CASE (g=7°) 


Table XII. 


£2300" 

% increase 

A/B 0 no tilt ον tilted (decrease) 
2.00 26,542.0psi. 23,006.0psi. 13.30 
2.33 26,666.0 22,662.0 15.00 
2.40 23,247.0 20.056.0 13.30 
2.66 26 ,669.0 22,167.0 16.90 
2.80 23,499.0 19,298.0 17.90 
3.00 26 ,680.0 22,213.0 16.75 
3.20 23 ,421 0 19,198,0 18.00 
3.50 26,629.0 210 19.00 
4.00 26,775.0 20,817.0 22.20 
^ E % increase 

A/B y no tilt y tilted (decrease) 
2.00 17,590.0psi. 14,990.0psi. 15.30 
2.33 17,657.0 14,962.0 15.30 
2.40 18,344,0 15,427.0 15.90 
2.66 17,699.0 14,697.0 17.20 
2.80 18,366.0 15,053.0 18.00 
3.00 19,407.0 16,102.0 17.00 
3.20 18,374.0 15,040.0 18.20 
3.50 19,492.0 15,656.0 19.70 
4.00 19,544.0 15,170.0 22.40 


Lateral and tilt stresses vs. 


aspect ratio (a/b), 7° tilt angle. 


Clamped approximation. 





— 


L:700" 

E f 5 | % increase 
A/B y no tilt y tilted (decrease) 
2.00 14,634.0psi. 12,373.0psi. 15.60 
2.93 14,667.0 12.459 0 15.60 
2.40 14,713.0 12,346.0 16.10 
2.66 17,699.0 14,679.0 17.20 
2.80 14,814.0 12,112.0 18.30 
34 00 14,843.0 11,918.0 19.75 
3.20 14,700.0 11 6070 20.60 
3.50 14,913.0 11,664.0 22.50 
4.00 14,884.0 Ll. S100 25.10 


Table XII. Lateral and tilt stresses vs. 
aspect ratio (a/b), 7° tilt angle. 
Clamped approximation. 
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Figure XXII. Stress magnification vs. 
aspect ratio (a/b), 7? tilt angle. 


Clamped approximation. 
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EDGE E 
L:300' 

E انك‎ i 9 ' ΜΟΝ | 
A/B (£23 % increase (g=5”) % increase 
2.00 25,023.5psi. 5.75 24,013.2psi. 9.45 
2.33 24,946.0 6.40 23,902.5 10.75 
2.40 21,877.0 5.90 20.965.5 9.30 
2.06 24,735.0 7.30 23,449.0 12.10 
2.80 24,762.0 7.20 23,485.0 11.985 
3.00 21,695.0 7.66 20,494.0 12.80 
3.20 21,607.9 7.74 20,401.1 12.70 
3.50 24,528.3 7.90 23,041.6 13.50 
4.00 24,216.8 9.58 22,514.0 15.90 
L:500' 

O y me m ول‎ tilted f 
A/B (@#=3`) % increase (ὅ55”) % increase 
2.00 16,335.0psi. 6.60 15,666.0psi. 10.90 
2.33 16,499.0 6.60 15,729,0 10.90 
2.40 17,091.4 6.85 16,258.0 11.40 
2.66 16,402.0 7.30 15,530 12.20 
2.80 16,943.2 7.75 15,996.5 12.90 
3.00 17,948.0 7432 17,042.0 12.20 
3.20 16,942.0 7.78 15,989.0 13.00 
3.50 17,844.29 8.48 16,748.0 14.10 
4.00 17,665.0 9.63 16,415.6 16.00 


Table XIII, Lateral and tilt stresses vs. 
aspect ratio (a/b), 3° and 5° tilt angle. 
Clamped approximation. 
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L:700' 
“y tilted ۱ ὃν tilted f 

A/B )2-3 ( % increase (Y=5”) % increase 
2.00 13,663.0psi. 7.10 13,016.0psi. 11585 
2433 13,663.0 Tato 13,016.0 11.25 
2.40 13,697.4 6.95 13,028.3 11.5 
2.66 13,669.4 7.50 12,943.0 12.30 
2.80 13,653.0 7.84 12,881.0 13.05 
3.00 13,586.8 8.43 12,751. 1 14.05 
3.20 13,402.0 8.35 12,539.0 14.70 
3.50 13,474.4 9.67 12,517.0 16.10 
4.00 13 ,434.4 9.75 12,471.0 16.20 


Table XIII. 
aspect ratio (a/b), 3° and 5° tilt angle. 


Lateral and tilt stresses vs. 


Clamped approximation. 
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Figure XXIII. Stress magnification vs. 
aspect ratio (a/b), 3° and 5° tilt angle. 


Champed approximation. 
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V. DISCUSSION OF RESULTS 


Upon examination of the results of the computer study 
for all fifty-four cases, it was apparent that the effect 
of web tilt, while present to some degree at any point in 
the plate, had real significance only for the stress in the 
"y" direction (o, ) - That this should be so seemed intuitive 
at the outset, since the vector of the applied tilt moment 
parallelled the "x" direction of the plate panel. The re- 
sults confirmed this, and hence the values of lateral and 
tilt stresses in the "x" direction (0,) have not been pre- 
sented in the results. 

Another intuitive conclusion, that the major effect, 
or tilt stress variation, would occur along the central 
plate axis x=a/2, normal to the vector of the applied moment, 
was also borne out by the results. Figures VII and IX show 
the total stress (ον) variation along the y=0 axis for both 
support approximations. The maximum stress magnification 
occurs at the plate center in each case. This is also the 
point of maximum stress due to lateral load in each case, 
arising from the assumption that the plate was simply sup- 
ported along its transverse edges. Consideration of Figures 
VI and VIII, total stress (ον,) patterns along the x=a/2 
axis, shows that the magnitude of the tilt stress is a maxi- 
mum at its point of application, and diminishes across the 
plate axis. Since the additional stress is largest along 
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the x=a/2 axis as well, the discussion will be limited to 
this effect. 

Figure X and XV are presentations of the stress 
arising from web tilt alone, shown along the axis x=a/2 for 
several cases and angles of tilt. Each pair of plots com- 
pares the effect produced under similar geometry for the 
two support approximations. Consideration of these figures, 
as well as of the accompanying tables, IV and V, reveals that 
the actual magnitude of the tilt stress is very nearly in- 
dependent of the support conditions assumed for the plate 
panel. This says, in effect, that the shear pressure at 
the piate edge does not vary measurably with the edge sup- 
port conditions. This fact is not obvious from an examina- 
tion of equations (34), (13), (14), (21), and (22), the 
analytical expressions for the shear pressure coefficient. 

The chief conclusion to be drawn from the state- 
ments above lies in the fact that the major factor in stress 
amplificotion due to web tilt, the shear pressure, can be 
considered a constant under the most widely varying support 
conditions normally made for plate analysis. 

Zt is worth noting at this point, that the actual 
values of tilt stress obtained in this analysis must be con- 
sidered low when extrapolated to an actual ship plating 
problem. To demonstrate the reasoning behind this statement, 


consider again the analogy to a simple beam under a 
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distributed load (q). For a single span, analogous to the 
single panel examined in this thesis, the shear force at 
either edge is of the order of "qa/2," whereas for a double 
span, analogous to a continuous panel system, the shear 
force at the central support is of the order of "ga." Thus, 
in a real case, the tilt stresses might be expected to ex- 
ceed those found in this analysis. 

Before any discussion of the results presented in 
Figures XVI through XXIII, it is pertinent to discuss the 
non-dimensional parameters against which these results are 
presented. In most work which has been done on plate panels 
and systems, the fundamental identifying characteristic of 
the panel is the aspect ratio (a/b). In this work, it was 
apparent from the initial assumption of the magnitude of 
the tilt moment arm, (c SIN f), that the value of "s" would 
have an important bearing on the results. Additionally, 
since the plate stresses are inversely proportional to the 
square of the plate thickness, (t), the non-dimensional 
depth ratío (c/t) seemed pertinent. Analysis of the input 
data, all selected by standard design techniques, revealed 
that a rough but generally consistent relationship exists 
between these two ratios (Figure V), that is, as the plate 
aspect ratio increases, an increase in the depth ratio is 
to be expected as well. Since both of these ratios are 


readily obtainable in an engineering situation, (€ can be 
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taken to be equal to the full depth of the stiffener with 
very little sacrifice of accuracy), the results shown in 
Figures XVI through XXIII have been presented against both 
of these parameters. Any other non-dimensional relationship 
for further analysis may be obtained by reference to Table 
III. The results shown are piotted for each assumed ship 
length, showing, in the lower figure of each group, the 
corresponding design stress levels, The difference in 
allowable stress magnitude arises from the design technique 
employed [10] which permits a higher level of plate bending 
stress in a shorter ship, where the girder stress would 
presumably be smaller. 

Figures XVI through XIX show results at the panel 
center for the simple support approximation, while figures 
XX through XXIII show results at the center of the longi- 
tudinal edge (x=a/2, y=b/2) for the clamped approximation. 
In both situations the magnification has been presented at 
the point of maximum plate stress due to the lateral load 
alone. The fact that the stress megnification shown for 
the clamped approximation is ín all cases the greater arises 
from two facts: 

l. The stress at the panel center under a simple 
support approximation is greater than the stress which ob- 


tains at the plate edge under the clamped approximation, and, 


T 


(o c-elaw pee d axe Ais 0 9 (mm e hn Aa 
A A < 
— -— A As 
¡mr rm e AS an 
w — p x A u 
TE mu w - — al m a R 
mulis dd UL Ὅν αὶ θα 
وه‎ 54d — cap com | 
ce——— ^ c col -— o — — 
a u :. τ. καν α ειἰκαακὲ © ow fom 
aut —— u^ o Pa w , —  * | 


— 
TIE =-  ὟὝ Leni! TE eom 2 


- | D ÉL ÉL AAA ww e A 
P -— = 2 ۷۳۲ ۱ ار رسيب ۳۷۷۷۶ عبت سس‎ 

















ΙΕ. ۱ "- | — Go oc o قافتا‎ 7 
11 E A اس واس ا‎ 4 td) de Jof eg | 
- Iie AS M. ll ια 
| jai b. í. A NP 
= = = ια | Us o8 MN. 

(ro ue - 
' — | S aaa Ü  —a— w -- . 
- — iit 3i kk اا‎ 


æ - =- á- οτι... 


οτε» 


2. the magnitude of the tilt stress remains es- 
sentially constant for varying support conditions, as haS 
been pointed out, and decreases nearly linearly as the 
distance from the plate edge increases (Figure X through 
xv). 

Of particular interest in the results is the fact 
that a fairly regular relationship is implied between the 
value of stress magnification and the plate aspect ratio 
(Figures XVIII, XIx, XXII, and XXIII). Thus, while still 
not defining the absolute value of the tilt effect, it may 
be concluded that the effect of web tilt increases in im- 
portance for a particular plate panel as the aspect ratio 
of the panel itself increases. 

The question which this thesis seeks to answer asks 
not only the magnitude of the tilt effect, but, more funda- 
mentally, whether the effect is large enough to be of con- 
cern in surface ship construction. In order to properly 
define an answer to this question, it is necessary to con- 
sider first the nature of the design method used for scantling 
selection on surface ships. To a very simplified first 
approximation, the plating thickness of a particular panel 
is selected to withstand the lateral load or to resist in- 
stability failure in compressive loading, whichever is the 


more limiting. The stiffener scantlings are selected 
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principally to withstand the girder bending of the ship, that 
is, the chief considerations are usually the area of material 
in the stiffener and its section modulus. The choice between 
a number of stiffeners which all satisfy the required values 
of section modulus is generally based upon considerations of 
weiaht or such characteristics as instability of the web 
itself. In any event, the ratio of the depth of the stif- 
fener to the plating thickness (c/t) has been merely a de- 
rived value, never a design consideration. 

From Figures X through XV, it can be seen that the 
maximum value of tilt stress, for the conditions examined, 
is measured in thousands of pounds per square inch, and even 
this, as has already been discussed, is to be considered a 
low value. Stresses of this magnitude axe certainly not to 
be ignored in any design consideration. 

As a general conclusion then, it may be stated that 
initial web tilt can indeed have an important effect upon 
plate bending stresses. This is especiaily true upon con- 
Sideration of the improvements presently being made in 
structural design. As design techniques improve, effects 
which have been heretofore hidden within the “factor of 
safety” will become increasingly important. Construction 
methods which allow the depth of a stiffening member to 
be increased and its orientation to the plate changed for 


such reasons as support for an item within the hull or for 
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a more favorable connection to another member are to be con- 
sidered suspect. In such cases, an estimate of the effect 
οξ the chance might be obtained by reference to a plot such 


as is shown in Figure XXII or X. 
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VI. CONCLUSIONS 


Recounting the conclusions developed in the pre- 
vious section, they are 

1۰ The value of tilt-induced bending stress is es- 
sentially independent of the support conditions of the 
plate. 

2. The effect of an initial web tilt increases in 
a regular manner as both the panel aspect ratio (a/b) and 
depth ratio (¢/t) increase. 

3. The ratio of stiffener depth to plate thick- 
ness (c'/t), where c' can be taken as very nearly equal to 
£, should be considered a design parameter, and its effect 
accounted for in any case of an intentionally skew stif- 
fening member. 

4. The magnitude of tilt-induced stress, while not 
catastrophic, is large enough to merit further considera- 
tion in any design technique which seeks to optimize the 


structure, i.e., reduce the factor of safety in the design. 
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VII. RECOMMENDATIONS 


Having concluded that the effect of initial web tilt 
can be of importance in structural design under the assump- 
tions made in this thesis, the logical recommendations fol- 
low that 

l. the investigation be extended, and, 

2. the assumptions in the present investigation be 
themselves investigated, and experimentally verified. 

in line with these general statements, specific 
recommendations are put forth as follows: 

1. An analysis of ship design procedures should be 
made to determine what general relationship between the aspect 
ratio (a/b) and the depth ratio (c/t) exists. If no good 
relationship can be found, the use of the depth ratio (¢/t) 
as a basic design parameter for cases of skew stiffeners 
should be further developed. 

2. The present investigation should be extended to 
stronger materials and higher lateral loads to verify the 
consistancy of the data presented herein. 

3. An analysis of the actual deflection behavior of 
an initially tilted stiffening member under lateral load 
should be undertaken to determine the effect of such a de- 


flection upon the actual angle of tilt along the member. 
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Appendix A 
Detailed Solution of Leteral Load Case, 


Clamped Approximation 


In this case, the boundary conditions are, as men- 


tioned before 





»= 0, 9-220 at x = 0, a (5) 
Ox 
4 | 3 
f JO d 
d. LIE =F + (2-v) - | aty=+b/2 (6) 
AX T X oy 
2% م‎ at y = + b/2 (9) 
dy | 一 


The fundamental equation of the deflection surface which 


must be solved is 


2 afa acu 
ἘΝ + 2 —25. 3 + CU = q/D (26) 
IX ox oy dy 
or 
V * = q/D (2a) 


Assuming a general solution form of 

eom mq ۳ م‎ (48) 
where 5 is the displacement of a simply supported strip 
parallel to the x-axis, then, 
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(49) and V “e, = 0 (50) 


co 


V ^u, E 


For the case of a uniform load, as can be shown in [5], 


— 3 3 3 
uM τ S (x —2a x” +a x) (51) 


and, in order to solve for Ways let 


_ MIX 
دلا‎ — X. SIN E (3) 


mzl,3,... 


where, 
ga” 7 
y = == a cos Y + y ΚΕ som YY ب+‎ c sim €f 
m D Lm a m a a m a 
a 
+p BY cosy BY | (4) 
m a a 


This form of solution was developed by Nádai [7] and Levy 
[8], and has been presented by Timoshenko [5]. 

To evaluate the constant terms in Ya’ it can be seen 
that the deflection surface under a purely lateral load will 


by syrmetric, that is 
a(x,b) = w(x,-b) (52) 


hence the equation for the deflection surface must be an 
even function of the variable "y" by considerations of sym- 


metry. The terms "y COSH(y)" and "SINH(y)" are in themselves 
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odd functions and hence the coefficients of these terms 


must vanish, i.o., 


C =D = 0 (53) 


Thus, the solution can now be expressed as 


uo Dy + قا‎ 
q í 4 3, 3 aa” τ mmy 
= San τε -22 x +a X | + D 1 (A „COSH - 


m=1,3,... 


(54) 


this can also be expressed as a single series, viz., 
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da c | | 
oe , (—— +A com τε DOY sim PY jgpy تلا‎ 
be m^ 7? m a m a a a 


m-1,3,... 
(55) 

By inspection, it can be seen that this equation 
satisfies boundary condition (5), since both the sine term 
and its second derivative vanish at x = 0, a. Further, since 
symmetry has been considered, the remaining boundary condi- 
tions need only be satisfied at y = + b/2 to complete the 
solution. 

Differentiating the deflection equation above, and 


substituting into boundary condition (6), gives 
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4 0 [ 7 
Lt As L (v-1)SINH À) A COSH Mm | 
m T 
tay) -vicosa ap] 
+ B o (l+v=a, A )SINH A,+A, (1-V)COSH A, (56) 
where 
_ mob 
^n 7 28 m 
and 
a a=. اتلك‎ (16) 


ra EIS 


as before. 


Differentiating again, and substituting into boundary 


condition (9) gives 
0 = A SINH À tB, (SINH LR COSH ha) (57) 


Solving these two expressions simultaneously for the 


coefficients Am and Ba’ gives 


4 r “On (SINH | high mC OSH hee) E 


| (58) 
E nq? L a à +a SINH À COSH A +2SINH? A J 
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u: A | qa. SINH ^m | (uo) 
m 5 5 L — 
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If the constant term 5 S is factored out of these two 
m T 
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expressions, the final solution for this case can then be 


written in the following form. 








4 ο, 
4qa x = | 3 
w = 1 一 | 1*A, COSH AY +p BY stn RC ¡SIN Rx 
5 m m a m a a a 
T D =l ,3;,... 
(12) 
where, 
TT a SINH À FA y SOSH Aa un; 
m 5 2 u 
mmm FO SINE À „COSH 2 SINH 5 
and, 
ca SINH Amn 
Β E COEM Mg (22) 
3 } 1 
Cma Opp SERE À „COSH Am ^SINH An 
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Appendix B 


Detailed Solution of Tilt Load Case 


As stated in the discussion in the body of the paper, 
the boundary conditions applicable to this case are 


32m 5 
الل‎ = 0, ME g 0 


x 








^ at x = 0, a (5) 
4 r 3 = 
El ٩ ۵ ے‎ D | LE. (2-v) 9 e | at y = + b/2 (6) 
4 i 2 3 
ax dy dx “dy 
2 8 
-D | 2 +. 2 ع = أ‎ (x) at y = + b/2 (10) 
oy OX 
2 . 2 
-D ۱ M + ν 5 7 | = f, (x) at y = -b/2 (11) 
Ay Ox 
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Assuming a general deflection form, as in the previous cases, 
of 


- nx 
Sm” Yn SIN s (3) 
m=] 

which satisfies boundary condition (5) and where 

y =a SIN W + p cosu W% + c sm I 

m m a m a m a a 

+p EY cosy SY (4) 
m a a 


This is the same form of solution technique as was used in 


the solution of phase one, the difference in the expressions 
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for Y, arises from the fact that there is no lateral loading 
to be considered in this case. 


Considering first a symmetrical loading form, in 








which 
M, = A, (60) 
y= +b/2 y= -b/2 
and assuming that 
£, (x) = f, (x) 5 En SIN T (61) 
m=1 


now, from consideration of symeretry, the odd functions present 


in V must pe eliminated, and nence 


A = D. = 0 (62) 
Then 
0 
= ` (n cosn BDO + c BY sım 2 ) sry UA (63) 
Z m a [a a a a 


m=] 


Having made use of symmetry to eliminate the odd functions, 
this deflection surface need satisfy the remaining boundary 
conditions at y = +b/2 only. 

Differentiating this deflection equation, and sub- 


stituting into (6) gives 


B = €, (to) (64) 
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where 


(1+v-0a A, )SINH hy 7 (lv) A COSH * 


- (1-v)SINH à +a COSH A (25) 
m m ra 
and, as before 
Q RUD „ti mr 
T 2n (15) and dae, (16) 


The deflection surface now can be reduced to 


a=} cip cosg D + WY com BE | sin BE (64) 
G WL a a a J a 
m=1 


gives 
* a < 1 
= = ο- |---- — E 
Cm D Gr R; (65) 
where 
Ry = | o (1-7) +2 COSH Aja PA gy (17) SINH Am (27) 


and hence, for a symmetrical moment loading, 
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(66) 


Considering now a parallel case in which the loading 


is completely unsymmetric, i.e., 
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y=b/2 y= -b/2 (67) 
and 
— EF MIX, 
f, (x) £, (x) و‎ Br SIN 2 (68) 
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from symmetry, it can be seen that, in this case, the coef- 


ficients of the even functions of "y" must vanish, i.e., 


B = Co = Q (69) 
and hence 
K 
Es ra Ty. τιν ny mix 
T L A, SINH š το, 2 COSH h ) SIN ich (70) 
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Proceedings as before, differentiating and substi- 


tuting into (5) 


A, T D. (v) (71) 
whore 
E (v*l-a, A, )COSH X -(l-v)A SINE A — 
Y SINH X +(l-v)COSH A 
m m 10 
and then 
i r ^ 5 
حا‎ oD fy sium 2. EE وون‎ E stn BE (72) 
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Differentiating again, and substituting into (10) gives 





E - 
=. القع ع‎ 二 
D D ) R- (73) 
where 
R, = Ly (l-v) +2 | SINH A 4A. (1-V)COSH À, (28) 
Thus, for an unsymmetrical loading, 
2 اس‎ n τ nay ^ IT 
«ow oh! LY grey BE ν BY cosy HA ' cry BE 
2 ALES lU EFE a aR a ل‎ a 
T D am 2 2 
msl 
(74) 


In the completely general case both of these solu- 


tions may be combined, as 
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2 A 
ys- -2 | L: L cosy Y y MIX guy BUY 
TD سه‎ m LR, a aRy 
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+ LL sim E + BY cosy EX syy ME (75) 
R a aR, ae 3 a 


since, under an arbitrary moment, all of the deflection 


functions, both odd and even, must be considered acting on 


the plate. 
To specialize this general solution for the case in 


which the moment is applied only along one edge of the plate, 


the given moment functions may be divided into both sym- 


metrical and unsymmetrical components, as 
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My - Ἂν = 5 L £4 (x) *£, (x) | (76) 
y=b/2 y= -n/2 

M = Mo >| £, (x) -£, (x) | (77) 
y=b/2 y= -b/2 





and, representing these separate functions as series, 


My =) Et SIN Pe (78) 
m=1 

M" = ۲۳ E" SIN u (79) 
m=1 


Now, if the applied moment is to be distributed only along 


the edge y = + b/2, then, 


£, (x) = Q (80) 
and 
E! =p” - ᾱ- (81) 
τὰ m 2 n 


thus the solution becomes 
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Appendix C 


Evaluation of the Second Order Shear Term 


In this analysis, the magnitude of the tilt moment 
acting on the plate is a direct function of the shear pres- 
sure or reaction which the plate exerts on the stiffener 
under a lateral load. There is, however, another shear 
pressure present, that arising from the tilt moment itself, 
which could add or detract from the initial value, depending 
upon the direction of application of the tilt moment, and 
so change the overall effect on the plate. 

In order to gain some insight into the magnitude of 
this increment of shear, it is helpful to consider a simple 
beam under lateral load, and compare the form of the result 


to those obtained in the body of the paper. 
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Figure XXIV Simple Beam Analogy 





Since, in the example pictured, the value of the shear force 
used in computing the tilt moment would have to be corrected 


by the amount of shear it itself created in the plate, it 
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can be seen that a corrected moment would be, when carried 


out through several steps, 


_ Whe W ς, E as Rau 3 
T, = > SIN Y + 5 (c SIN £) F 21 (£ SIN 2( + ..٠ 


)82( 


By inspection, it can be seen that the relationship between 


terms is 


c SIN Y (83) 


L 


Considering a numerical example with values consistent with 


those used in this investigation, assume: 


L = 72 inches 
€ = 7 inches 
g = 7 degrees 


the ratio of the first two terms is then, 


112.22) = 0.015 


or, the error in neglecting this tem is of the order of 
1.5%. The following terms decrease rapidiy and can be 
neglected as well. 

Considering the case in point, it has already been 
shown that the shear pressure (reaction) at the plate edge 
is 


— م‎ sIN EZ (84) 
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where A is a series constant developed from equation (32) 
and arises from the lateral load alone. From this has been 


derived the tilt loading factor, E ， 





4ga 
E = FP, € SIN Y (85) 
m rm m 


The shear loading on the plate from this, i.e., the moment 
load, alone can be obtained by differentiating (24) and sub- 


stituting into the expression, 
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This results in 
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or, simplified by grouping the series constant into a tem 
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E بت‎ 时 En Ro, SIN 5 (87) 
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hence, the order of magnitude of the increment of tilt 


moment is 
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Again, this is related to E, by the ratio 


T ¢ SIN 2 
2a (89) 
as was true for the beam analogy. This correction is then 
of second order in (e SIN g)? as well, and may be neglected 


with little error. 
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Appendix D 


Computer Programs 


1. Program for the stresses in a plate subject to 


lateral load. This program is presented for both of the sup- 


port approximations. In addition this program calculates 


the value of the coefficients of the shear pressure term (E). 


* 








> NM... see = Eee ee 


FORTRAN 
STRESSES IN A LATERALLY LOADED FLAT PLATE 


INTERNAL FUNCTIONS 
SIHHP(X) - (ExPF(X)-EXPF(-X))/2.0 


COSHF (X) = (8XPP(X)+EXPF(-xX))/2.0 


— تا ی را‎ oco. MIROR 





o | Q 





FOR THE SIMPLE SUPPORT APPROXIMATION 
AAMF (X,Y) -(EN*(1.0--EN) *SIBHF (X) -EH* (1.0-EN) *COSHF (x} 

1 -Y*(2.0*cosur (x) »x*sIBUF(X)))/((3.0*EN) * (1. 0-EB) 

2 *siNiP(x)*cOsHF(X)-(1.0-EN)**2*X-2.0*Y*(COSHF (X) **2)) 
BBMF(X,Y)-(EN*(1.0-EN)*SINHP(X) Y*COSHP(x))/((3.O0-*EN) 

1 *(1.0-EN) *SIMHF (X) *COSHF (Xx) - ) 1 , 0-511 ( **2%*x-2,08y8 


2 (cOSHF (X) **2)) 


E .. ae Ee Rae DF sr um 








FOR THE CLAMPED APPROXIMATION 
AAMP (X,Y)sY*(-SINHF(X)-x*COSHF (X))/(x*v«v*sINnr (X) 
J *COSHF (X)2.0* (SENHF (X) **2)) 
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BBMF(X,Y) Y*SINUE (X)/(X*Y+Y*SINHF (X) *COSHF (X)+2.0 
* (SINHF (X) **2)) 

READ COUNTER FOR NUMBER OF DATA RUNS 

LITEM=0 

READ 920, INDEX 

LITEM=LITEM+1 

READ INPUT DATA 

READ 900,SHIP,YIELD,A,B,H,S,Q,EN 

PRINT 921,SHIP,YIELD,A,B,H,S 

PI-3.141592654 

FRONT-(24.0*Q*A**2)/(H**2*PI**3) 

ONE-(12.0*S*(1.0-EN**2)IPI)/(H**3*A) 

SECTION TO CALCULATE STRESSES ALONG THE X-A/2 AXIS 

11=0 

DO 250 L=1,11 

I=L-1 

Y=FLOATF (I) *B/20.0 

SET UP TO COMPUTE TEN SERIES VALUES 

AMDA=P1I *Y/A 

SIGMAX=0.0 

SIGMAY=0.C 

DO 200 M=1,19,2 

FM=FLOATF (M) 

FIRST=(SINF (FM*PI/2.0))/(FM**3) 


AX=FM*ONE 
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200 


250 
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Cl-zFM*PX*B/(2.0*A) 
PL=AAMF (C1,A%) 
P2=BBMF (C1,AX 
XTERM=1.0+P1*(1.0-EN) *COSHF (PM*AMDA) -2, 0 *EN*P2 
*COSHD (FM*AMDA) -+P2* (1. 0-EN) *FM*AMDA*SINHF (FM*AMDA) 
UPl-XTEDM FIRST "FRONT 
SIGMAX=SIGMAX+UP 1 
VALUES OF STRESS IN THE Y DIRECTION 
YTERM=EN+ ( (EN-1.0) *P1-2.,0*P2) *COSHF (FM*AMDA ) 
+ (EN-1.0) *P2*FM*AMDAFSINHF (FM*AMDA) 
UP2=YTERM*FIRST*FRONT 
SIGMAY=SIGMAY+UP2 
CONTINUE 
SIGYX (II) =SIGMAX 
SIGYY (11) =SIGMAY 


CONTINUE 


SECTION TO CALCULATE STRESSES ALONG THE Y O AXIS 
11=0 

DO 350 L=1,11 

11-1 

II=II+1 

X=((PLOATF (1) *A)/20.0)+A/2.0 

VALUES COMPUTED AT TEN EQUAL STATIONS ON Y AXIS 


SIGMAX=0 .0 
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300 


350 
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SICMAY-0.0 
pO 300 M=1,19,2 

ZM=FLOATF (M) 

2LEAD=(SINP (2M*PI*X/A))/(2M**3) 
227 *ONE 

Z1-ZM*PI*B/(2.0*A) 
ZllzAAMF(Zl,AZ) 

212z-BBMF (21,AZ 
ZTERMz-Z1.04211*(J1.0-EN)-2.0*EN*z12 
UP3-ZTERM*ZLEAD*FRONT 
SIGMAX=SIGMAX+UP 3 

STRESS IN THE Y DIRECTION 
YITEM=EN+Z211*(EN-1.0)-2.0*Z12 
UP4=YITEM*2LEAD*FRONT 
SIGMAY=SIGMAY+UP4 

CONTINUE 

SIGXX (II) =SIGMAX 

SIGXY (11) =SICMAY 


CONTINUE 


PRINT HEADINGS 


PRINT 923 


PRINT 901 
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400 


450 
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DO 400 N=1,11 

L=L+1 

PRINT 912,L,SIGYX(N) ,SIGYY(N) 
CONTINUE 

PRINT 907 

J=9 

DO 450 N=1,11 

J=J+1 

PRINT 913,J,SIGXX(N),SICXY(N) 


CONTINUE 


CALCULATION OF COEFFICIENTS FOR TILT MOMENT 


PRINT 915 

DO 500 N=1,19,2 

FN=FLOATF (N) 

C3=FN*PI*B/(2.0*A) 

AN=FN*ONE 

P3=AAMF (C3,AN) 

P4=BBMF (C3,AN) 

EM=4 .0*Q*A (( (EN-1.0) *P3+ (EN+1.0) *P4) *SINHF (C3) 
+ (EN-1.0) *P4*C3 *COSHF(C3)) /(PI**2*FN**2) 


PRINT 909,N,EM 


CHECK FOR OTHER CASES 
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900 
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IF(INDEX-LITEM) 15,15,5 
CALL EXIT 


FORMAT STATEMENTS 


FORMAT (F5.1,F10.1,3F3.5,F3.2,2F8.5) 
FORMAT (35X, 29HSTRESSES ALONG THE X=A/2 AXIS//) 
FORMAT (1H1,35X,27HSTRESSES ALONG THE Y=0 AXIS//) 
FORMAT (10X, 2HM=,I2,5X,5HE(M)=,E20.9//) 

FORMAT (10X,2HY-,13,4HB/20,5X,9HSIGMA (X)z,E20.9,5X, 

OHSIGMA (Y) =,E20.9//) 

FORMAT (10X,2HX=,13,4HB/20,5X , 9HSIGMA (X) =,E20.9,5X, 

9HSIGMA (Y) z,E20.9//) 

FORMAT (1H1,20X,21HTABLE OF COEFFICIENTS///) 

FORMAT (12) 

FORMAT (151, ,7HLENGTH=,F5.1,5X,6HYIELD=,Flo.1/30X, 
2HA=,F8.5,5X,2HB=,F8.5,5X,2HH=,F8.5,5X,2H1=, 
F8.2////) 

FORMAT (30X,41HSTRESSES IN A LATERALLY LOADED FLAT 

PLATE//) 


END 


2. This section of coding is designed to receive 


the output from the first section already presented, examine 


the tilt induced stresses for all angles up to seven degrees, 


and compare the results. 
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FORTRAN 


TILTED STIFFENER 


DIMENSION EM(25),SIGYX(50),SIGYY(50),SIGXX(50), 
SIGXY (50) ,‚SIGTYX (50) ,SIGTYY (50) ,SIGTXX (50), 
SIGTXY (50) , TSIGYX (50) , TSIGYY (50) , TSIGXX (50), 


TSIGXY (50) ,SLFYX (50) , SLFYY (50) ,SLFXX (50) ,SLFXY (50) 


INTERNAL FUNCTIONS 


SXNHF (X) =(EXPF (X) -EXPF (-X))/2.0 

COSHF (X) =(EXPF (X) +EXPF (-X))/2.0 

RHOF (X, Y) =((1.0+EN- (Y*X) ) *SINHF (X) - (1.0-EN) *x* 
*COSHF (X) ) / ((1.0-EN) *SINHF (X) +Y*COSHF (X) ) 

GAMF (X, Y) =( (1.0+EM-=(Y*X) ) *COSHF (X) - (1. 0-EN) *x 
*#SINHF (X)) / (Y*SINHF (Xx) + (1.0-EN) *COSHF (X) ) 

ADENF (X, Y) =(RHOF (X, Y) * (1.0-EN)+2.0) *COSHF (X) +X 
*(1.0-EN) *SINHF (X) 

BDENF (X, Y) =(GAMF(X, Y) *(1.0-EN)+2.0) *SINHF (X)+X 
*(1.0-EN) *COSHF (X) 

READ 910,INDEX 

ICOUNT=0 


LCOUNT=LCOUNT+1 
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READ 900,A,B,EN,EPS,S, (EM(I) ,I=1,20) ,(SIGYX(I), 
I=1,24) , (SIGYY (LI) ,I=1,24), (SIGXX(I) ,I=1,24), 
(51οΧΥ (τ) 151,24) 

OUT-3.0/ (H**2) 

PI=3.141592654 

AMDA::(PI*B)/(2.0*A) 


ALPH-(12,.0*S*(1.0-EN**2)*PI)/(A*H**3) 


CYCLE TILT ANGLE FROM ONE TO SEVEN DEGREES 
DO 800 K=1,7 


THETA=(FLOATF (K) *PI)/180.0 


STRESSES ALONG THE A/2 AXIS, 21 POINTS 
Y=-B/2.0 
DELTAY=0.0 

DO 250 I=1,21 
SIGTX=0.0 
SIGTY=0.0 
Y=Y+DELTAY 

DO 200 M=1,19,2 
FM=FLOATF (M) 
X1=FM*AMDA 
AX=FM*ALPH 
P1=RHOF (X1,AX) 


P2=GAMF (X1,AX) 
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R1L=ADENF (X1,AX) 

R2=BDENF (X1,AX) 

XTERM=((P1*(EN-1,0)+2.0*EN/R1+FM*PI*Y* (EN-1,0)/ 
(A*R2))*cosHF (FM*PI*Y/A)-(FM*PI*Y*(EN-1.0)/ 
(A*R1)+(P2*(EN~1.0)+2.0*EN) /R2) *SINHF (EM*PI*Y/A) 

FRONT=EM (M) *SINF (FM*PI/2.0) *EPS*SINF (THETA) 

YTERM-((Pl*(1.0-EN)-2.0)/nl-FM*PI*Y*(1.0-EN)/ 
(A*R2) ) *COSHF (FM*PI*Y/A) + (FM*PI*Y(1.0-EN) / 
(A*R1)+(P2*(1.0-EN)+2.0)/R2) *SINHF (FM*PI*Y/A ) 

UPX=OUT*FRONT*XTERM 

UPY=OUT*FRONT*YVERM 

SIGTX=SIGTX+UPX 

SIGTY=SIGTY+UPY 

CONTINUE 

SIGTYX (1) =SIGTX 

SIGTYY (1) =SIGTY 

DELTAY=3/20.0 


CONTINUE 


STRESSES ALONG THE Y AXIS, 21 POINTS 
x=0.0 

DELTAX=0.,0 

DO 350 Iz1,21 

SIGTX-0.0 


SIGTY-0.0 
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X=X+DELTAX 

DO 300 J=1,19,2 

FM=FLOATF (J) 

X1=FM*AMDA 

AX=FM*ALPH 

P1=RHOF (X1,AX) 

P2=GAMP (X1,AX) 

R1=ADENF (x1,AX) 

FPRONT=EM (3) *SINF (FM*PI*x/A) *EPS*SINF (THETA) 
XTERM=(P1*(EN-1.0)+2.0*EN) /R1 
YTERM-(Pl*(1.0-EN)42.0)/R1l 
UPX=0UT*FRONT*XTEPM 
UPY=OUT*FRONT*YTERM 
SIGTX=SIGTX+UPX 
SIGTY=SIGTY+UPY 

CONTINUE 

SIGTXA (1) =S1GTX 
SIGTXY (1) =SIGTY 

DELTAX=A/20.0 


CONTINUE 


SUM STRESSES AND GET RATIOS 
ALONG A/2.0 AXIS 
DO 450 N=1,21 


TSIGYX (N) =SIGTYX (N) +SIGYX (N) 
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SLEYX (N) =TSIGYX (N) /SIGYX (N) 

TSIGYY (N) =SIGTYY (N) +SIGYY (N) 

SLFYY (N) =TSIGYY(N)/SIGYY (N) 

ALONG Y AXIS 

TSIGXX (N) SSIGTXX (N) +SIGXX (N) 

SLFXX (N) =TXIGXX (N) /SICXX (N) 

TSIGXY (N) =SICTXY (N) +SIGXY (N) 

SLFXY (N) =TSIGXY (N) /SIGXY (N) 

CONTINUE 

PRINT OUT VALUES FOR THIS ANGLE OF TILT 
PRINT 901,K,A,B,H,EPS,EN,S 

PRINT 902 

DO 600 N=1,21 

L=-11 N 

PRINT 903,1, TSIGYX (N), SLFYX(N), TSIGYY (N), SLFYY(N) 
CONTINUE 

PRINT 904 

DO 700 N=1,21 

L=N-1 

PRINT 903 ,L, TSIGXX (N), SLFXX (N), TSICXY (N), SLFXY(N) 
CONTINUE 

CONTINUE 

IF (INDEX-LCOUNT) 20,20,30 


CALL EXIT 
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FORMAT STATEMENTS 

FORMAT (6F10.5/(4E20.9)) 

FORMAT (1H1,50X,15HANGLE OF TILT=,13,7HDEGREES//5X, 
2HA=,F10.5,5X, 2HH=,F10.5,5X,4HEPS=,F10.5,3HNU=, 
F6.4,2HI=,F10.5//) 

FORMAT (45X,29HSTRESSES ALONG THE x=A/2 AXIS,/7X, 
5H20V/B,15X, GSHSIGMA (X) , 18X,6HFACTOR,18x, 
SHSIGMA (Y) ,18X, 6HFACTOR) 

FORMAT (7X,14,7X,E20.9,5X,F20.9,5X,E20.9,5X,E20.9) 

FORMAT (//,45X,27HSTRESSES ALONG THE Y=0 AXIS,/7X, 
5H20X/A,15X,8HSIGMA(X),18X,6HFACTOR,18X, 
SHSIGMA (Y),18X,6HFACTOR) 

FORMAT (I2) 


END 
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STRESSES IN A LATERALLY LCABED FLAT POSTE 


- NC RCTATICN AT ECCE 


STRESSES ALCNG THE X=A/2 AXIS 


Y: 000 20 SIGMA(X)= ۰521761447۲ 04 SIGMA(Y)= 11 
Y= 1B/2C SIGMA(X)= .511CC2€60E 04 SIGNA(Y)- .121101363E 
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Y= 8B/20 SIGMA(X)= -.1€5810406€ 04 SIGMA (Y)= -:105317420E 


Y= 98/20 SIGMALX)= - , 24 7826432۲ 04 5 1 6۳۵ ) ۷ ( < - , 1 6 62 9 2 1- 


Ys 10B72C — SIGMA(X)= 一 1999092351 4 5 ۲6۳۸ ) ۷ ( > - و‎ 2 342 1 782 


es 


۱1 
y 
ERII 
i 
is 
۱ 
I 
(9 1 
7 | 


۱ ۶ 





10A/2C 


114/20 


124/2€ 


134/20 


L4A/2C 


15A/2C 


16^/2C 


17A/2C 


18A/20 


19A/20 


202/0 


SIGMA(X)= 


SIGMA(X)= 


SICMA(X)= 


SICMA(X) = 


SIGMA(X) = < 


SIGMA{X)= 


SIGMA(X)= 


SIGMA{X)= 


SIGMA(X)= 


SIGMA(X)= 


SICMA(X)= 


9۱۳۳ Ες REG 


«221761447E 


۰ ۶215/0764 ۲ 


«5193659406 


ع 3 2 5244064 ۰ 


" م‎ 53 7EÉ 


ح 210 192.2 0 24 ۰ 


۰۶7432321 ۲ 


۰62 6040 4 ۴ 


e 6 2 9۱32 4 0 


«517292417E 


oe 2ETOEEEB2E- 


THE Y=0 AXIS 


04 


04 


04 


04 


04 


04 


04 


C4 


C4 


04 


02 


SIGMA(Y)= 


SIGMA(Y)= 


SIGMA(Y)= 


SIGMA (Y)= 


SIGMA(Y)= 


SIGMA (Y)= 


SIGMA (Y)= 


SIGMA(Y)= 


SIGMA(Y)= 


SIGMA (Y)= 


SIGMA(Y)= 


۰1 24697 


.124843946€ 


«124690130€ 


.124721640€ 


¿124035531E 


.122588511E 


»118138900E 


.«108760871€ 


¿890050985E 


»547697902€ 


05 


05 


05 


05 


05 


05 


05 


05 


04 


04 


.212254321E-02 


spomi) © ۷ 1 ٩۱ 4 

1 ۱ ۶ ۱ ' 
۱۱ 9 I > 449—200 ١ 

۱ ۷ ۶ ۲۲ = 41105 
۱۱۲ I 5 
Il 1 ie 
1,۶ - u "1 
RII 2! 
۱۱ ۷ 1 
MEX! - ‘to œ 

۱ - LT 

y 





“= 1 
ν- 3 
ν- 5 
M- 7 
ν- 9 
ν-11 
v=13 
M=15 
v=17 


TABUE OF 


E (M) = 


E(M)= 


EiM)= 


ELM ) = 


E{M)= 


EfM)= 


ElM)= 


ElM)= 


E(M) = 


E{M)= 


UPEPFPFTCTENINS 


SIRO IE 
.TEB3AITACE 
ου σον E 
2207 794638 
۰ 1. ٩5 ۶ E 
.8556099CIE 
|.61300122CE 
.4€0549772E 
.35861522CE 


.281124212E 


e 


02 


C2 


C2 


02 


Cl 


Cl 


Cl 


01 


ο] 








ANGLE OF TILT = 3DEGREES 


A= 96.00000 B= 30.00000 H= 20000 EDO T.T8000NU- .30001- 37.80000 


STRESSES ALONG THE XsA/2 AXIS 


SOVAS - SIGMA(X) 1 FACTOR SIGMA (Y) FACTOR 
-10 - 54569 1565E 04 .990440734E 00 -. 233007757E 05 .994833574E 00 
-9 -.339790225€E 04 .976895936E 00 -.164272875E 05 .987848215E 00 
-8 -, 154994 916E 04 .934771955E 00 -.1025052B4E 05 .973298475E 00 
7 ۰863957494۴ 2 -.1 75361 32 89] 01 -.478471432E 04 .930011950E 00 
-6 .150957689FE 04 .112088820E 01 -.342388302E 02 .723111652E-01 
-5 .271847565E 04 .107495530E 01 .399941228E 04 .114912793E 01 
-4 .371245004E 0% .106173538E 01 .731562503E 04 .108926998E 01 
-3 .449111424E 04 .105687354E 01 .991406195E 04 ۰107374972] 1 
-2 .505424075E 04 . 105576411۴ 01 ۰1179451626 5 . 106919523E 01 
-] .540168218E 04 ‚105707515E 01 .129568528E 05 ۰106991800۲ 1 
0 . 553335555۲ 04 . 106051446۲ 1 ۰ 1340100326۲ 5 ۰107467746] 1 
1 ۰544922531] 04 . 106637905E 01 . 131269442E 05 ۰108 3963376۶ 01 
2 ۰ 514930934۲7 4 ۰1 075622676 01 ۰1 21346938۲ 5 ۰110003298] 1 
3 .463367376E 04 .109042143E 01 .104243129E 05 .112901282E 01 
4 .390245002E 04 ‚11160 7408E 01 .799592294E 04 .119056387E 01 
5 .295584816E 04 -116881850E 01 .484970871E 04 .139343867E 01 
6 ۰179423012۴ 04 Ἢ .1 23224837۲ Οἱ . 985969618E 03 -. 208233200E 1 
7 .418199353E 03 - ,848837979E O1 -. 359476022۴ 04 .598718823E 00 
8 -.117117641E 04 -706334688E 00 -.889120676€ 04 .844229452E 00 
9 -.297233570E 04 .854545668E 00 -.148996277E 05 .895983018E 00 
10 -.498431228E 04 .90466230^E 00 -.216078898E 05 .922555298E 00 

| STRESSES ALDNG THE | 

20X/A SIGMA(X) * FACTOR SIGMA (Y) FACTOR 
0 .260068879E-02 ‚100C00000E οἱ .212254319E-02 . 100000000E οἱ 
1 .533412829E 04 .103116294^4E 01 ۰5367775071] 04 ۰103665739] 1 
2 ۰657353401۲ O04 . 104385841۲ 01 ۰9291 06720] 04 ۰1063880346 1 
3 .659577943E 04 .105362639E 01 .114307920E 05 .105100225E Ol 
4 .609859720E 04 .106187560E 01 .124979515E 05 .105790315E 01 
5 .583501298E 04 .106452742E 01 .130374659E 05 ۰106351451 ۶ O01 
6 ۰ 5591 14501۲ 04 . 106548425E 01 .132473923E 05 .106803205E 01 
7 . 557645574۲ 04 .10633843?2E 01. .133587228E 05 .107108301E Ol 
8 .551657617E 04 .106225496E 01 .133817885E 05 .107320352E Ol 
9 .557471350E 04 .106042336€E 01 .134112028E 05 .107423734E 01 
10 .553335555E 04 .106051446E O1 .134010036E 05 .107467746E 01 
jet .557471350E 04 .106042336E 01 .134112028E 05 .107423734E 01 
12 .551657617E 04 .106225496E 01 .133B17BB5E 05 .107320352€ 01 
13 .557645574E 04 . 106338432۴ 01 .133587228E 05 .107108301E 01 
14 .559114501E 04 .106548425E 01 .132473923E 05 . 106803207E 01 
15 .583501898E 04 .106452742E 01 .130374661E 05 . 106351451E Ol 
16 .609859720E 04 . 106187560201 .124979517E 05 .105790318E O1 
17 .659577951E 04 . 105362640۳ 01 ۰114307922۲ 5 ۰1051002266 01 
18 .657353416E 04 .104385842E 01 . 9291 06750E 04 .104388037E 1 
19 .533412859E 04% .103116299€ 01 .567775108E 04 .103665745E 01 
20 .293595087E-02 .112891279E 01 .250889767E-02 .118202431E O1 
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SIGMA (X) 


-.542187534E 04 
-.334443640E 04 
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e1/6654115E 03 
.161789544E 04 
.284458969E 04 
۰ 385606661 ۶ 4 
۰465 1907 35۲ 04 
.523185171E 04 
.559572443E 04 
.574342236E 04 
.267489862t 0% 
.539017066E 04 
. 488931358E 04 
«417247601E 04 
°. 323988922E 04 


.209193211E 04 


.729210921E 03 
.847217187E 03 
.263573494E 04 
. 463484276E 04 


SIGMA(X) 


H 


.544137917E 04 ` 


.675728783E 04 
.681912847E 04 
.633502677E 04 


.607033737E 04 


۰581976630] 04 
.579760015E 04 
.573167637E 04 
.578604981E 04 
.574342236E 04 


.578604981E 04 


.573167644E 04 
.579760015€E 04 
.581976630E 04 
.607033737E 04 
.633502677E 04 
۰6819128626 ۶۸ 
۰675728805۲ 04 
۰5441 37970۲ 04 
.315900505E-02 
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EPS= 


FACTOR 
. 984080844E 
. 96 T32 4509C 
.B91374908€ 


-. 35856277 76 


.120131668E 
.112482402E 


 .110280874^E 


.109471225E 


,109286468E 


.109504802E 
. 11007157596 
.111054190E 
. 112593543E 
.11505 7998E 
.11932996 7E 
.128113564^E 
. 15532974 9E 


-.148011209E 


‚510955378E 
۰ 7 < 7 77304۲ 
. 941 2329 1 1 


FACTOR 
.100000000E 
.105189607E 
.107303798E 
. 108930472E 
.110304224E 
.110745838E 
.110905181E 
.110555474€ 
.110367399E 
.110062379E 
.110077553E 
-110062379E 
.11036 7401 
.110555474E 
.110905181€ 
.110745838E 
.110304224E 
.108930476E 
.107303802E 
.105189615€ 
.121468015E 
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STRESSES ALONG THE X=A/2 AXIS 


Y=0 AXIS 


SDEGREES 


. 30001 = 


SIGMA (Y) 
۰. 2382 02689۳ 05 
۰ 16292 84361 5 
۰ 100634336۲ 5 
, 45451 5:2 75۲ 04 
.258002613E 3 
.434472512E 04 
.17145096T7E Οὐ 
.103671003E 05 
.123023544E 05 
.135201845E 05 
.140205514E 05 
.138034398E 05 
.128688565E 05 
.112168279E 05 
.884741858E 04 
.576073483E 04 
.195696835E 04 


-.256350677E 04 


.779973T60E 04 


-.137488119E 05 


۰204010837] 05 


SIGMA (Y) 
. 2122543198 -2 


. ۰581132680] 04 


۰955091022۲ 04 
.117998444E 05 
.129530666E 05 
.135554885E 05 
.138088094E 05 
.139485620E 05 
.139890701E 05 
.140278203E 05 
.140205514E 05 
.140278203E 05 
.139890701E 05 
.139485620E 05 
.138088096E 05 
.135554887E 05 
.129530668E 05 
.117998446E 05 
.95509107^E 04 
.581132740E 04 
.276594434E-02 
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.991396286€ 
.979 63485 
.955533624E 
.88344B012€ 
.544892140E 
۰21 248 3446 7 
۰11 + +5 
۰112281 
۰ 1۳ 19 
۰11 164358 ع‎ 
۰11 2+36 ( ۲ 
۰ 1 139 ۲ 
.116658623E 
٠121484573 
.131734846E 
.165519852E 
۰ +13 3047۲ 
.498 2725 72 
. 740593 07 
.826779179E 
.871030353E 


FACTOR 
.100000000E 
.106104603E 
.107307452E 
.1084934T1E 
.109642688E 
.110577152E 
.111329467E 
.111837544€ 
.11219067TE 
.112362840E 
۰ 11 24361 ۴ 
.112362840E 
۰ 2 121 906 77 
۰1 11 83 ۲ ۲ 
۰1 1 1 69 ع‎ 
۰ 1 1۳57 ۶1 
ه‎ 1 0642690 
۰ 10849347 3 
۰1 0 ۲ 
۰ 1 0610461 5 
۰ 1 30 312745 





A 


= 96.00000 


20Y/B 


! 
E | | | ١ | (١ | >» 
m MW DUO WOOO 


O Y OU YOU T Q í = O 


= 


20X/A 


PP Ó) fN ©‏ + 0ن پد ο‏ و 


8 


30. 00000 


SIGMA (X) 
.538694188€ 04 
.329113357E 04 
.140625495E 04 
.266637225E 03 
.172588367E 04 
.297031909Ε 04 
.399924520E 04 
.481221005E 04 
.540892102E 04 
.578917503E 04 
.595284842E 04 
.589988381E 04 
.563029751E 04 
.514417373E 04 
.444167856E 04 
.352306403E 04 


.238872629E 04 


.103927401E 04 
.524245940E 03 
.230016064E 04 
.428643901E 04 


SIGMA(X) 
.260068879E-02 
.554830298E 04 
„694048 114E 04 
. 704179645E 04 
.657073528E 04 
.630493812E 04 
.604769044E 04 
.601807006E 04 
.594612069E 04 
.599674165E 04 
.595284842E 04 
.599674165E 04 
.594612069E 04 
.601807006E 04 
.604769036E 04 
.630493812E 04 
.657073535E 04 
.704179659E 04 
.694048144E 04 
.554830365E 04 
.338137906E-02 


H 


.50000 


ANGLE OF TILT = 


EPS= 


FACTOR 


-977740355E 


.946199961€ 
.848110199E 


.128149990E 
,1174540TO0E 


.114375683€ 


.113243558E 
. 112985212 
hls 290508 E 
. 114091382 
. 11545 ۴ 
. 11 ۲ 
۰1 2 10۶ 7 
ολο ουσ. 
.139311027E 
E تاد رت‎ 


- . 21 ) 946 1 5 ۲ 


. 2 191۳۷ TSE 
۰ 6612 
, / 7 7 996 95 


STRESSES ALONG 


FACTOR 
.100C00000F 
«ΠΟ. Ὲ 
«110212856ε 
.112487428€ 
.114408333€ 
.115025841E 
. 11524864 غ7‎ 
. 114759655 
.114496673€ 
.114070166E 
٠. 114509 1382 
. غ114070166‎ 
. 1 144966 غ73‎ 
.114759655E 
. 115248647 
.1150258^1E 
.114408335E 
.112487429E 
. 110212861 
غ107256607.‎ 
.130018594€ 


THE 


01 
01 
01 
Ol 
01 


7.78000NU= 


STRESSES ALONG THE X=A/2 AXIS 


Y=0 AXIS 


TOEGREES 


«300015 


SIGMA (Y) 


-.231400065E 05 


.161588097E 05 
.987590911E 04 


. 949352862E 03 
. 468898498E 04 
۰ 91121 7۲785۲ 04 
۰1081875 71 6 5 
۰128086438 5 
۰ 1 4021 791۲ 05 
.146382097E 05 
.144778723E 05 
.136007803E 05 
.120069259E 05 
© 969631746E 04 
. 666898265E 04 
۰ 2 92500708] 04 


-.1535398 18E 04 


.671159692E 04 
.126015056E 05 
.191979576E 05 


SIGMA(Y) 
.212254319E-02 
.594449557E 04 
.980996080E Οὐ 
.121677713E 05 
.13405793BE 05 
.140719315E 05 
.143685147E 05 
.145366022E 05 
.145944994E 05 
.146425573E 05 
.146382097E 05 
.146425573E 05 
.145944994E 05 
.145366022E 05 
.143685147E 05 
.140719319E 05 
.134067941E 05 
.121577718E 05 
.980996154E 04 
.594449639E 04 
.302220710E-02 


37.80000 


FACTOR 
.987969480E 
.971703410E 
.937822938E 
.837026067€ 
.116021329E 
۰.1 4 ۲ 
.120787380E 
.117173339E 
.116112784€ 
.116281088E 
.117389373E 
.119551688€ 
.123293651E 
۰1 3004+ 6 
.144374643€ 
۰1 91 ۵1601 
.617750868E 
.298437599E 
.637273200E 
.7571786378E 
.819662526E 


FACTOR 
.100000000E 
+. 108536030E 
.110217964E 
.111876369E 
.113483313€ 
.114789970E 
.115841926E 
.116552366E 
.117046149E 
ell7286884E 
.117389373E 
.11728688^E 


.117046149E 


.116552366E 
.115841926E 
.114789973E 
«1134833 17E 
«111876372E 
۰ 11021 97 
۰ 1 08 6045 
۰ 1+2 3861 
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